Phase transition to bundles of flexible supramolecular polymers 
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We report Monte Carlo simulations of the self-assembly of supramolecular polymers based on a 
model of patchy particles. We find a first-order phase transition, characterized by hysteresis and 
nucleation, toward a solid bundle of polymers, of length much greater than the average gas phase 
length. We argue that the bundling transition is the supramolecular equivalent of the sublimation 
transition, that results from a weak chain-chain interaction. We provide a qualitative equation of 
state that gives physical insight beyond the specific values of the parameters used in our simulations. 



Self-assembly is an active field of research, driven 
by the desire to design new materials. Understand- 
ing the rules of self-assembly has been defined as the 
challenge for chemistry of this century []]]. The large 
molecules involved in self-assembly spontaneously orga- 
nize because of highly specific interactions like hydropho- 
bic/hydrophilic, hydrogen bonding, Coulomb interaction 
and 7T — 7r-stacking. Modeling these self-assembly pro- 
cesses using coarse-grained models, such as the patchy 
particle[2], has the potential to yield theoretical insight. 
Up to now, patchy-particle models have mostly been used 
to describe the self-assembly of functionalized colloids for 
photonic crystals [3K and to study the formation of self- 
assembled clusters[4|. 

Here, we focus on linear supramolecular self-organi- 
zation, such as the reversible aggregation of aromatic 
molecules and discotic proteins into supramolecular poly- 
mers, which has been the subject of several experimen- 
tal studies^, @, 0]. In this context, Sciortino et al.[2] 
have shown that the first-order Wertheim Thermody- 
namic Perturbation Theory (WTPT) of associating liq- 
uids accurately predicts the chain-length distribution by 
comparing the WTPT to simulations of simple square- 
well patchy particles. One step beyond linear polymer- 
ization is the association of supramolecular chains into 
bundles 0,0,0], a problem that neither the WTPT, nor 
the square-well patchy particle, can deal with, since it 
requires chain-chain interactions. 

In this Letter we present Monte Carlo simulations of 
a type of patchy particle that, by decreasing tempera- 
ture or increasing density, first polymerizes into chains 
and subsequently undergoes a phase transition toward 
bundles of these chains. We interpret this bundling as a 
sublimation transition from a polymer gas to a solid bun- 
dle. This sublimation competes with polymerization and 
gives rise to non-trivial phase behavior. We propose a 
simple thermodynamic model to describe the transition. 

Our model is a coarse-grained representation of a disc- 
like, aromatic molecule and consists of hard spheres 
of diameter a, dressed by two opposing patches. The 
orientationally-dependent patch potential allows not only 
for chain formation, but also exhibits a weak chain-chain 



Figure 1: a) Geometry of the patches. 0% is the angle between 
the direction pi of patch i and the interparticle vector ry = 
fj — fj. 0j is the angle between the direction pj of patch j 



and 



6ij is the angle between pj and —pt. b) 



Patch potential for three values of 0j for w = 0.4. Notice that 
the minimum at —1 occurs for 6i = 0j = 6ij = 0, and that 
it decreases rapidly for increasing angles. The particles are 
depicted at half their size for clarity. 



interaction. The patch potential between patch i directed 
along the vector p, and patch j with direction vector pj , 
illustrated in Fig. Q] is given by a Lennard-Jones poten- 
tial of the interparticle distance r = \fj\ modulated by 
three directional components 
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where e is the maximum energy of interaction, and w pe- 
nalizes non-perfect alignment. We truncate V(r) of Eq.Q] 
at r c — 2.3a and shift and rescale it to have V(r c ) = 
and y(2 1 / 6 cr) = — e. The first two directional compo- 
nents favor minimization of the angle 9i between the 
patch direction pi of patch i and the interparticle vec- 
tor fij, and of the angle 9j between p"j and fji = —fij. 
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The third component minimizes the angle 9ij between pj 
and — pi, favoring parallel alignment of the patches. This 
potential has three advantages compared to a square- 
well patch potential^. First, it allows multiple bonds 
without increasing the energy per interaction site in dis- 
crete steps. This feature introduces a slight interaction 
between molecules in neighboring chains. Second, the de- 
sired parallel alignment of neighboring patches prevents 
branching (and therefore network formation) of the poly- 
mers. Finally, it allows us to tune the flexibility of a 
supramolecular chain. One can show that w is a measure 
of the chain flexibility, by calculating the bending rigidity 
k, and hence the persistence length l p of a chain 
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where (3 = (fceT) -1 , L is the chain length, and R(s) 
is the chain radius at s. For example, the oligo(p- 
phenylenevinylene)-derivative OPV-4 in dodecane has 
l p = 150nm at 300K, which, for a molecule separation of 
a = 0.35nm and a bonding energy of 56kJ/mol[6] yields 
w w 0.2. 

We study the thermodynamic equilibrium of a sys- 
tem of M patchy particles in a simulation box of volume 

V with periodic boundary conditions at temperature T 
by Monte Carlo simulations. The system is equilibrated 
by performing moves and rotations of single molecules, 
and of whole chains. In addition, we perform reptation 
moves. We note that we can simulate a limited tem- 
perature range (fc^T > 0.04e) as at lower temperatures 
the probability of removing a molecule from a chain by 
a simple Monte Carlo move vanishes. We choose to con- 
sider molecules bonded if they interact with an energy 

V < — 0.3e. Contrary to the square-well patch poten- 
tial [2] , the potential of EqQ] makes the average energy 
per bond -Ebond temperature dependent. We find that 
equipartition (i.e. -Ebond = — e + iksT) applies up to 
k B T < 0.07e where E hond 0.79e. 

Supramolecular polymerization is reversible and there 
exists an equilibrium density p(N) of chains with length 
N, with N = 1,2,..., oo. The total density is given by 
p = J2n=i Np(N). The average chain length N, and the 
aggregation fraction T] are defined as [t| 
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The fraction of sites that are not bonded X, is given by (3] 



X = 



T,n=i 2N p( n ) 



(3) 



where the factors 2 appear because each molecule has two 
possible bonding sites. In the ideal free-association (IFA) 
model[3| the polymers form an ideal gas, and each bond 



Figure 2: The aggregation fraction r\ (Eq. D as a func- 
tion of temperature for several values of the flexibility w for 
M — 1330, p — 1.13- 10 -3 . r) = 1/2 defines the polymerization 
temperature T*, indicated by an asterisk on the correspond- 
ing curves and on the horizontal axis. At w = 0.4 a transition 
from a gas of chains, to a solid bundle at k B T = 0.057e is in- 
dicated by the dotted line. The symbols denote simulation re- 
sults, the solid line is the WTPT. For comparison the dashed 
line shows a fit of the IFA-model to the curve of w = 0.3, with 
v = 1.4 • 10~ 4 cr 3 and G = 0.91e. 



lowers the energy by a discrete amount e. This results in 
p(N) ~ p{^) N , and gives 
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N = - + -y/l + 4pvexp(l3G), 



(4) 



where G < e is an effective free energy per bond 
and v is a bonding volume. Both parameters are not 
known a priori, and are usually fitted to the average 
chain lengths, determined by e.g. circular dichroism 
measurements Q. The WTPT includes the spatial ex- 
tension of the molecules [13] neglected in the IFA-model, 
by calculating a reference hard sphere free energy and 
adding the attractive contribution of the pair potential. 
WTPT also assumes that no rings can be formed and 
that only one interaction per attractive site is possible. 
The average chain length predicted by the WTPT is for- 
mally equal to the IFA-result of Eq. [4] if 



4tt J 5rep (r) (exp(-/3K tt (r)) - 1)^ r'dr 
J exp(/?G) 



(5) 



where A involves a single site-site interaction, and is re- 
lated to the second virial coefficient [2], (.) Wl)W2 denotes 
an average over all orientations uj\ and u>2 of the two 
molecules, flrepM is the pair correlation function of the 
hard-sphere part of the potential, and V^ttM is the at- 
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tractive part. At low densities and temperatures the IFA- 
model and the WTPT are equivalent, and in that case, 
for the square- well patchy particle of Ref. [2f , G is simply 
the well depth and v can be calculated analytically. 

In Figure 2 we compare the aggregation fraction 77 as a 
function of temperature T for several values of the flexi- 
bility w to the prediction of the WTPT with g Ye p{ r ) = 1, 
appropriate at low densities, and A from Eq. \E\ numer- 
ically calculated. The remarkable agreement with sim- 
ulations shows that the WTPT also holds for smoothly 
varying potentials on a hard sphere. For comparison we 
also show a fit of the IFA-model, that deviates at higher 
temperatures due to the temperature dependence of the 
association energy. We find (not shown) that this de- 
viation reduces with decreasing w, and is negligible for 
w = 0.1. The polymerization temperature T*, defined as 
the temperature where half of the molecules in the system 
has aggregated, i.e. where Tj(T*) = s, increases with w. 
This rise in T* is due to an increase in available bonding 
volume v, or, equivalently, because a transition from an 
unbound to a bound state costs less entropy for a more 
flexible chain. For w = 0.4 and UbT < 0.058e the ag- 
gregation fraction suddenly jumps to rj ~ 1. The chains 
have bundled, with a concomitant increase of the aver- 
age chain length and a strong depletion of the gas density. 
Such a sudden increase of aggregation, not accounted for 
by polymerization theory, has recently been observed for 
the OPV-4 molecule Moreover, the transition is rem- 
iniscent of one of the assembly pathways suggested for 
the formation of amyloid fibrils [ll| . 

In Fig. 3 we examine the bundling transition in more 
detail. Figures 3a and 3b show the dramatic differ- 
ence between a bundled and a polymer gas configuration, 
while the temperatures differ by only fc^AT = O.Ole. 
In the bundle the individual chains remain identifiable, 
with no bridging connections between chains. Visual 
inspection of configurations in the process of bundling, 
shows that when three chains come together, they remain 
bonded and suddenly grow in length, suggesting a nucle- 
ation mechanism. In Figure 3c the average chain length 
N is shown as a function of temperature for w = 0.4, 
where the bundling transition is visible as a sharp jump. 
The dashed line represents systems heated up from a 
bundled configuration at ksT = 0.057e. This hystere- 
sis together with the nucleation mechanism, is evidence 
for a first order phase transition. We identify the ac- 
tual bundling temperature T& at the high end of the 
hysteresis- loop, at fcsTf, m 0.065e. Nevertheless, only 
at ksT = 0.057e the critical nucleus is small enough to 
appear spontaneously during the duration of one simu- 
lation. In the inset of Fig. 3c we show JV as a function 
of T for several values of w. Although a stiffer chain 
loses less entropy upon bundling than a flexible one, Tb 
increases with flexibility w, similar to T*. The transi- 
tion is thus not driven by entropy, but by the interaction 
energy of neighboring chains (Eq. [TJ that, at large 9, 




T{e/P) 

Figure 3: a) Top panels: representative configuration for 
w = 0.4 at two very close temperatures. M, p as in Fig. 2. 
a) ksT = 0.057e, where a solid bundle has nucleated, b) 
ksT — 0.058e where the system is still a polymer gas. c) 
Average chain length W as a function of temperature for 
w — 0.4. The solid line shows the jump of the chain length 
for ksT ~ 0.057e on cooling. The dashed line shows the hys- 
teresis of the transition while raising the temperature from an 
equilibrated system initially at fesT = 0.057e. Inset: Same 
curve for several values of w. No hysteresis loop is shown. 
Notice that the bundling temperature increases with w. 

increases with w. Assuming that the patches in neigh- 
boring chains are perfectly aligned, the smallest angle 0$ 
between molecules in neighboring chains is of the order 
9i ~ arctan(2~iJ) w 0.73rad or 42°. For w < 0.15 this 
lateral interaction is negligible (6i/2w w 2.4). Increasing 
w also increases the available lateral bonding volume, de- 
creasing the entropy loss upon bundling. 

We interpret the bundling as a sublimation transition, 
from a gas of polymers to a solid bundle. To derive an 
approximate equation of state, we equate the chemical 
potential /i so i of a bundle to that of a polymer gas /j, gas . 
Using the IFA-model we can derive 

Mga S = |§ = k B T In (^j - 2/T 1 m^V (6) 

where the first term on the right hand side is the chemical 
potential of an ideal gas of monomers with internal par- 
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Figure 4: The gas solid coexistence line (sublimation line) 
with G = 0.8e, Z\ = 1, e so iid = -1.3 for w = 0.4, v = 3e-3cr 3 , 
and p = 0.0013a" 3 fitted to k B T* = 0.065e resulting in k = 
2.1 • 10 3 ecr" 2 . The solid line is calculated from Eq. El the 
dotted line from [9] the dashed line from Eq. 1101 In the inset 
we show the same, but on a linear density scale. 



tition function Z\. For our rigid molecules, Z\ = \. As a 
first approximation, we model the bundle by an Einstein 
CTystal[12j with 



Msoi — e S oi + ksT 
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where e so i is the well depth, including the association 
energy G and the various chain-chain interactions of a 
molecule in the bundle. Furthermore, k is the spring con- 
stant of the Einstein crystal. We replace the de Broglie 
wavelength term with the bonding volume A 3 = v, as 
was done to derive the IFA-model. We measure e so i in 
a simulated bundle and fit k to the sublimation transi- 
tion of Fig. 3 (p = 3 • 10 3 cr~ 3 , k B T b = 0.065e). Solving 

Mgas = Msol 

yields the sublimation density Pg ub ((3) 



exp(-(3p, so] )Z 1 



v [Z 1 exp(/3G) - exp(-pp so] )Y 
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Equating the chemical potential of an ideal gas to p, SO ] of 
Eq. [7] results in the sublimation line for an ideal gas of 
monomers p^^iP) 
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exp(/3e so iid) 



(9) 



It is useful to also define the density at the polymerization 
temperature p* — p(T*), by combining Eqs. [2]and|4] 



V(T*) = - 



P 



v 1 exp(— (3G). (10) 



Comparison of p^ uh to p* allows us to estimate whether 
the sublimation transition is dominated by the bundling 
of polymers (Eq. [8]), or of monomers (Eq. [9]). 



In Figure 4 we compare the sublimation line of a poly- 
mer gas Pg ub to that of an ideal gas of monomers p™ b in 
coexistence with the same Einstein crystal, and the poly- 
merization line p* . At low densities pP uh ~ p™ lb , because 
the transition from gas to solid occurs at a higher tem- 
perature than polymerization, i.e. Xf, > T*. The subli- 
mation lines differ at higher densities, because p gas tends 
to the binding energy G, whereas the chemical potential 
of an ideal monomer gas tends to zero. At high densi- 
ties polymerization occurs at higher temperatures than 
bundling, i.e. Tt> < T* . Sublimation requires e so i < — G. 
As decreasing w lowers e so i , we expect a shift of the sub- 
limation line p su b to lower temperatures, increasing the 
polymer-dominated regime. When the lateral interac- 
tions vanish, i.e. below w = 0.15, the bundling transition 
will completely disappear. 

In summary, we have presented a model based on 
patchy particles, that describes supramolecular polymer- 
ization and displays a first order phase transition to bun- 
dles, due to weak chain-chain interactions. The chain-to- 
bundle transition can be seen as a sublimation transition 
from a polymer gas to a solid bundle, for which we have 
been able to give a qualitative equation of state. We 
have related the occurrence of the phase transition to 
the flexibility of the supramolecular polymer. Our simu- 
lations show that bundling leads to a sudden increase of 
the average length of the aggregates, as experimentally 
observed in different polymerizing systems, from OPV-4 
to amyloids. 

We believe that the patchy particle is a flexible, power- 
ful tool for efficient modeling of complex, self-assembling 
systems that can be adjusted at will by including more 
specific, e.g. chiral, chain-chain interactions. 
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